Fluctuations in a mesoscopic superconducting ring: resonant behavior of conductivity 

and specific heat in two mode critical regime 
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The critical fluctuations in a mesoscopic superconducting ring are studied within the Ginzburg- 
Landau approach. The nonlocal conductivity as well as the specific heat are calculated as functions of 
the magnetic flux $ through the ring. At $ = $o/2 two low-energy eigenstates become degenerate 
and near this point the behavior of fluctuations-dependent quantities change dramatically: both 
the zero Fourier component <to of the fluctuation conductivity and the specific heat become non- 
monotonic functions of | <& — <J>o /2| with rather special resonant structure. 
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I. INTRODUCTION 

Mesoscopic superconducting rings attract considerable 
interest in connection with their unusual properties in 
the vicinity of the superconduction transition. One of 
the most interesting features is an oscillatory behavior 
of the critical temperature as a function of magnetic 
flux through the ringi due to competition between spa- 
tial modes characterizing by different orbital momenta. 
The remarkable property of the superconducting rings is 
the fact that all eigenstates are well separated from each 
other and in some cases the dynamics of each state can 
be described independently. This provide a unique op- 
portunity to study the critical fluctuation contribution 
to all measurable quantities of the rings. With the devel- 
opment of the microfabrication technologies it becomes 
possible to create superconducting rings of the radius R 
comparable to the Ginzburg-Landau coherence length £o 
(see Refs. [HH and references therein). The fluctuations 
in such small rings can be described in the frames of 
simple analytical models, such as the zero-dimensional 
Ginzburg-Landau formalism (see, for example, Ref. 0). 
In Refs. y@ the persistent current due to the critical 
fluctuations was studied for different ratio between the 
rings' radius and the coherence length. The calculation 
of magnetization in the critical regime have been per- 
formed in Ref. 0. Note that previously the influence of 
the critical fluctuations on specific heat and magnetiza- 
tion has been studied in small superconducting granula 
both theoretically^ and experimentally^. 

At the same time one of the most natural way to study 
fluctuations near T c is the performance of the conduc- 
tivity measurements (the fluctuation correction to the 
conductivity has been described first in Ref. [Toh . For su- 
perconducting rings two types of experiments were pro- 
posed: (i) nonlocal paraconductivity measurements with 
probes attached to different points of the rin g n i 12 ; (ii) 
contactless measurements of dissipation in the array of 
rings subject to an electromotive force induced by a weak 
alternating magnetic fielcU^. The corresponding calcula- 



tions of the nonlocal conductivity within the Gaussian 
approximation were performed in Ref. [lil 

The separation of different eigenstates in mesoscopic 
ring with R > £ allows to describe analytically the con- 
tribution of the critical fluctuations to the nonlocal con- 
ductivity in magnetic field which produces the flux $ 
through the ring. It was shown that in the critical fluc- 
tuation regime the fluctuation contribution to the con- 
ductivity has a logarithmic singularity near the critical 
temperature with a prefactor proportional to the 
(^/^o) 2 - However this approach is valid only for mag- 
netic flux values which are not close to $o/2. Otherwise 
for $ sa $o/2 the separation of modes with different or- 
bital momenta breaks down: critical fluctuations near T c 
are produced by two interacting modes, which can not 
be considered separately anymore. 



In the present paper we suggest an analytical descrip- 
tion of critical fluctuations for the case $ « ^o/2 and 
calculate the corresponding fluctuation correction to the 
non-local conductivity of the ring. We will keep in mind 
the contactless realization of conductivity measurements 
in which only the zero Fourier component of the conduc- 
tivity plays the key role. Also we analyze the behavior 
of the specific heat value in two-mode critical regime. To 
analyze the situation we use time-dependent Ginzburg- 
Landau equation. This approach is reasonable for small 
dirty superconducting rings (see Ref. [l3h . 

The paper has the following structure. In SectionHIlwc 
describe in detail the behavior of the conductivity both 
inside and outside of the two-mode critical regime. In 
particular, we present the exact expression for the con- 
ductivity for the case <fi w 1/2, based on the analytical 
solution of nonlinear two-mode Ginzburg-Landau equa- 
tion. In Section Hill we calculate the specific heat taking 
the mode interaction into account. The results are sum- 
marized in Section ITVl 
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II. FLUCTUATION CONDUCTIVITY DUE TO 
CRITICAL FLUCTUATIONS 

To calculate the non-local conductivity we use the ap- 
proach which is similar to that of Ref. [lU The zero- 
frequency conductivity a(ip — i p') is given by the Kubo 
formula (see, for example, Ref. Il4l ) 



T 



o{ V - V ') = - \ (J(tp,0)Jfr',t))dt, 



(1) 



where the supercurrent is defined by the standard expres- 
sion J ((p, t) = (e/mR)Re [ip (id v - rp) ip*] and = $/$ - 
Performing the Fourier transform of the order param- 
eter wave function and the conductivity 



*l>(<p,t) 



4> n (t)e mv , a(cp-ip') = y"V fc e 



(2) 



one can obtain the following expression for the zero 
Fourier component CTq 



Gim) = y/2b/ a 2 T c V , where V = 2nRs is the volume of 
the ring and s is its cross-section. In the dirty limit a sim- 
ple estimate shows that Girp) <x (T c /E F )[(£ l) 3/4 /VV], 
where / is the electron mean free path. Ovcrwise when 
e* <C Giro) the fluctuations of one or two lowest modes 
become critical and nonlinear terms in Eq. ([5]) begin to 
play the key-role in the description of the fluctuation 
conductivity. Further we will consider these fluctuation 
regimes separately 



A. Gaussian regime 

In the Gaussian regime one can find from the lin- 
earized Eq. (O that ip n (t) = Vra(O) ex P(~ t/ T n), where 
t- 1 = {8T c0 /it) [e* + {£, /R) 2 n(n - 20)] . Also it is easy 
to find the exact expressions for correlators (|^„(0)| 2 ) 
and (l^nCO)! 4 ) which contribute to ao : 



<l^(0)| 2 ) 



4mf§T 



V [e* + {^/R) 2 n{n - 20)] ' 



(6) 



CT = 



T \2mR 



I ^]T(n-0) (m-0) 

q n.rn 

x(|0„(O)| 2 |V m (t)| 2 ), 



(3) 



while for nonzero components a^o 



oo 



+ k - 20) 



(4) 



The Fourier components of the order parameter wave 
function ipn{t) satisfy the non-linear time-dependent 
Ginzburg-Landau equation^ 



jd t ip n = a(T)ip n + b 



(5) 



k+m—p—n 



where the values a = (4m^) 1 [e+ (£,o/R) 2 (n — <j)) 2 ] and 
b are the parameters of the Ginzburg-Landau theory, e = 
(T - T c0 )/T c0 , £q~ 2 = 4maT c0 and 7 = 7ra/8. 

Since the critical temperature T c {4>) is a periodic func- 
tion with the period 1 due to the Little-Parks effect 
and T c (—<fi) = T c (cp) we can consider only the range 
< < 1/2 for the magnetic flux. In this range 
T c {4>) = T c q [l — (£0 / R) 2 2 ] and it is convenient to use 
the parameter e* = [T - T c (rf>)] /T c0 = e + {£, /R) 2 4> 2 
instead of e. 

Let us analyze the dependencies of the conductivity on 
the magnetic flux through the ring in different tempera- 
ture ranges. In case when when e* ^> Gi( ) the system is 
in the Gaussian regime and one can neglect the nonlin- 
earity in Eq. ([3]) . Here we have introduced the Ginzburg- 
Levanyuk number Gi( ) characterizing the width of the 
critical fluctuation region in zero-dimensional systems: 



<l^(0)| 4 ) 



2(4m^) 2 T 2 



V 2 [e* + (Zo/R) 2 n(n-2(f > )Y 



(7) 



From the expression ([3]) one can see that the singular 
part of cto comes only from terms with m = and m = 1. 
The corresponding expression for the fluctuation conduc- 
tivity in the Gaussian regime is rather complicated and 
therefore is given in the Appendix. Here we will focus 
only on two limiting cases. 

For e* <C (£o/R) 2 (l — 20) the most singular part in the 
temperature dependence of the conductivity comes from 
the mode with n — and has the form 



AG) 



{el 



Atts 2 



(8) 



For e* > (£ /i?) 2 (l - 20) in the vicinity of the half quan- 
tum flux through the ring two modes with n = and 
n = 1 contribute to the most singular part of cfq, which 
reads as 



AG) 



(£*) 



Airs 2 



2(1 - 20) 5 



(9) 



These expressions are valid for any 
0< 1/2. 



from the range < 



B. Critical regime with noninteracting modes 

Let us first consider critical fluctuations for magnetic 
flux values which are far from = 1/2. In this case 
only one spatial mode with m = is in a critical regime 
and thus different spatial modes do not interact with 
each other. Then the equation for %po(t) can be 
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solved analytically^ 3 -. Note that the dominant contribu- 
tion to the conductivity comes from the correlator with 
n = m = . Considering only this correlator in the 
expression ([3]) (the lowest mode approximation) one can 
reproduce the result of Ref. [IH with corrected numeri- 
cal coefficient (the corr elator <|^ (0)| 2 |Vo(*)| 2 ) should be 
calculated instead of the correlator (V'o(O)V'o(O) 2 ) which 
provide an additional factor (87r) 1 / 2 r- 2 (l /4) 0.38): 



JO) 



(0) 



r 3/2„2 



' 1 1 1 

• 3 In— 

G« (0) e* 



(10) 



Note that the lowest mode approximation, which has 
been used to obtain the expression (|10l) , is valid when the 
sum of all neglected terms in © is small compared with 

. Let us calculate all terms with m = and m = 1 in 
the expression ([3]) which can become singular near T c (0) 
for < 4> < 1/2. The corresponding expression for the 

( o 

conductivity <7q is given in the Appendix since it is 
rather cumbersome. 

Near = 1/2 for e* < (£ /.R) 2 (l - 20) the most sin- 
gular correction to the conductivity (|10[) is negative and 
can be written in the form 



(0) 



8tts 2 



1 



1 



« 2 0) (1 " 20) 



ml. 



(ii) 



From the expression (fTTj) one can see that for 

|0-1/2| < {R/Zof Gi (0) the value Acr<, 0) /a^ 0) exceeds 
1 which indicates that the expression (|10|) is not valid in 
this region. Note that in our case R > £o and the typical 
value for the Ginzburg-Levanyuk number is Gim\ ~ 10~ 3 
(for Al ringsii), so the lowest mode approximation breaks 
down only in a narrow region near the point = 1/2 
when both modes with m = and m = 1 are strongly 
fluctuating and interact with each other. The case when 
|0 — 1/2| < (i?/£o) 2 Gz(o) is considered in the next sub- 
section. 



C. Two-mode critical regime 

To obtain the expression for the critical fluctuations' 
contribution to paraconductivity at the magnetic flux 
w 1/2 one should consider the dynamics of two low- 
est interacting modes. Further we restrict ourselves to 
the most interesting case when the magnetic flux through 
the ring is exactly half quantum = 1/2. In this case all 
modes except the modes with n = 0, 1 can be neglected 
since they do not contribute to the singular part of (To. 
Then the expression for Co reads 

1 



T \2mR, 

x((l^o(o)| 2 -|Vi(o)| 2 ) / (\Mt)\ 2 ~\Mt)\ 2 )dt 



The system of nonlinear time-dependent Ginzburg- 
Landau equations for ipo (t) and 0i (t) is 



r ^o + 6 (|^o| 2 + 2|^i| 2 ) -00 = 0, 



4 ^o (12) 
7^1 + ^^2^1 + b W + 2 I^| 2 ) ^ = °> 

where e* = e + £q / AR 2 . Let us multiply the first equation 
of this system by ipQ, the second equation by Vi and 
then sum each of the obtained equations with its complex 
conjugated. The result is 



d~ t F + a*F + F 2 + 2FG = 0, 
<9rG + a*G + G 2 + 2.FG = 0, 



(13) 



where i = 16bt/na, a* = e* /Amt^b, F = \tpo\ 2 and G = 
\ipi\ 2 - To simplify these equations we introduce a new 
time variable r = (a*) -1 [l — exp(— a*t)] and new func- 



tions /(t) 



- p a,*t(r) 



F{t(r)) and g(r) 



Then the system Q13p transforms into the form 

d T f + f + 2fg = 0, 
d T g + g 2 + 2fg = 0. 



G(t(r)). 



(14) 



To obtain the exact solution of the system (fT4")l it is 
convenient to consider the auxiliary functions u = (/ + 
g)/2 and v = (f — g)/2. Note that Co depends only on 
v(t): 



(15) 



From the system (|14[) wc obtain the equations for new 
functions 

d T u = v 2 - 3u 2 , d T v = -2uv, (16) 
which can be reduced to the equation for the v function 



2vd 2 v - 5 (d T v) 2 + 4v 4 = 0. 
Then for the decaying solution of Eq. (|17l) we find 



d T v = -2(j,v 2 Vl + Xv, 



(17) 



(18) 



where A 



v 2 )/vq (here and after the index in- 



dicates the function value at r = and vq is assumed to 
be nonzero), = sign(«o)- 

Note that the explicit solution of the system (fT4| inte- 
grated over time has the form 



m 2 ^- go)2 for , = 1, 



f(t)dt 



ln- 



foa* 
9o 



In 



g(t)dt 



l/o - go\ 
fo 

l/o ~9o\ 



for 



for jx = 1; 



(19) 



In 2 ^^") 2 for n 
g a* 



4 




4V* 2(T-T(<t>))/Gi 



FIG. 1: (color online) The specific heat as a function of mag- 
netic flux through the ring at the temperature of the super- 
conducting transition for (R/^o) 2 Gi(o) = 10~ 2 . The exact 
solution (blue solid curve) near = 1/2 qualitatively dif- 
fers from the solution obtained taking only zero spatial mode 
into account (red dashed curve with circles) and is in sharp 
contrast with the predictions of the independent modes ap- 
proximation (black dashdot curve). 



An interesting feature of the solution (f!9| is that de- 
pending on the initial condition only one mode (/ or g) 
becomes slowly decaying due to the nonlinear modes' in- 
teraction. 

The expression (fT8|) allows us to calculate ctq since 



l/a* 



«(l/a*) 



v{r)dr 



(OrV) 



-dv. 



(20) 



To estimate the v(l/a*) value at the temperature of the 
superconducting transition (a* — > 0) we use the asymp- 
totic form of the equation ([15} solution, which gives us 
v (l/a*) ps /ia*/2. Then the singular part of the integral 
(IT%t has the form 



(21) 



The Gibbs free energy functional for two considered 
modes has the form 



F = V 



a* (/o+.9o) + - (/, 



9l 



4/oflb) 



(22) 



Performing the integration over f and g in the expres- 
sion ([15]) in the limit a* — s- we obtain 



cro( 



1/2) =A- 



7T 3 / 2 S 2 



J_l 1 

rv3 <=•* 



(23) 



FIG. 2: (color online) The dependencies of the specific heat 
on temperature for <f> — (red dashed curve) and 0=1/2 
(blue solid curve). The black dashdot curve corresponds to 
the prediction of the independent modes approximation. 



where A = n(V3-l) U/2arcth (^V^jJ ~ 1-42. Thus 
at = 1/2 the exact expression of Co is 1-42 times larger 
than the one obtained within the lowest mode approxi- 
mation (see the expression (JTOJ) ) . 

Note that for = 1/2 the nonzero components <7fc do 
not diverge at e* = 0. Indeed only the correlators in ((4]) 
that contain ipo or ipi functions can produce the singular- 
ity of crfc at e* = 0. But the prefactor of the correspond- 
ing correlator in the a± expression is zero in case 0=1/2 
and thus the <j\ component does not contain the singular 
part. To calculate other <Jk components we should take 
into account the following: (i) the correlators which con- 
tain noncritical modes can be represented as a product of 
two independent corelators; (ii) the corelators for critical 
modes can be calculated only for r = since near e* = 
the decay times t n = irR 2 /8T£ 2 n(n — 1) for noncritical 
modes are much shorter than for modes with n = 0,1; 
(iii) for noncritical modes (i/;„(0)^* (t)) = (\tp n (0)\ 2 ) t n . 
Note also that for 0=1/2 the correlators (Jo) = (go). 
Therefore one can easily make sure that does not con- 
tain a singular part at all since all contributions from 
critical modes cancel. Thus at = 1/2 only the <jq com- 
ponent has a singularity at e* = 0. 

It is interesting to analyze the conductivity behavior 
when the flux is slightly smaller than half-quantum <50 = 
1/2 — <C 1. Then the system of equation (TT4"]) has the 
form 



d T f + f 2 + 2fg = 0, 

A > Sa * ,2 

1 - a*r 



2/5 = 0, 



(24) 



where 5a* = 5<fi/ (2mR 2 b). Note that here the value a* 
which is used in the definition of functions / and g cor- 
responds to modified critical temperature and reads as 
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a* = (4m^)-i [e+(£ Q /R) 2 (l/2-66) 2 ]. In what fol- 
lows we will assume that 56 -C (i?/£o) 2 Gi(o)- 

For fo > go the dynamics of the slow-decaying mode /, 
which define the singular contribution to the conductiv- 
ity, is weakly affected by the differences in the character 
of the mode g decaying due to finite 5a*. At the same 
time for go > fo the situation changes. For the r values 
when g » Sa* /(l — a*r) the mode g is slowly decaying 
while at r ~ (a* + 6a*)' 1 which corresponds to the case 
g ~ 5a* /(l — a*r) the decaying of the (/-mode becomes 
exponential. Thus for go > fo one has to replace the 
cutoff I /a* in the integral flU) with l/(a* + 5a*). Then 
the expression (f2"3")l becomes modified as follows 

= ^fev (I) ( ln ? + ln ^w) ' 

(25) 

where <5e* = 2(£ /R) 2 56. 

It is interesting to compare the expressions (fTTj) . (f2"3")l 
and (|25l) . One can see that at magnetic fluxes close to 
6 = 1/2 the conductivity has rather special resonant 
structure. Indeed far from 6=1/2 the conductivity 
is well described with the expression (fTU|). When 6 in- 
creases the conductivity becomes smaller compared with 
the expression <|10[) since the mode with m = 1 begins 
to enter the critical regime. This picture is broken down 
when 1/2 — 6 « (i?/£o) 2 Gi(o), where the logarithmic 
singularity in conductivity (Tq(6) reveals a very sharp 
peak due to the nonlinear modes' interaction. The max- 
imum of the peak at 6 = 1/2 is larger than the one 
obtained in the lowest mode approximation (see the ex- 
pression (|I0[) ) . The width of the peak is decreasing while 
e* —> (see the expression (|2"5|l ). 

III. FLUCTUATION SPECIFIC HEAT 

Now let us turn to the calculation of the fluctuation 
specific heat AC (6, T). The explicit form of the Gibbs 
energy G[ip] for two interacting modes with n = 0, 1 was 
obtained in Ref. 0- Then the specific heat can be calcu- 
lated as AC = -Td\G. For \6- 1/2| > (R/tof Gi [a) 
the main contribution to the AC at e* —> comes from 
the mode with n = 0, which can be considered indepen- 
dently from other modes. Then 

AC( r = r cW) = fcS 5 l-( 1 -|^). W 

The corresponding dependence AC (6) at e* — is shown 
in Fig. 1 (red dashed curve with circles). Note that at 
6 = 1/2 the exact calculation gives the value AC = 
9 Gi-*[1 - (£o/2i?) 2 ], where 9 = (4/3) (1 + VSrj - ttt? 2 ), 

■q = (arcch 2)" 1 . 

Moreover it is interesting to analyze the dependence of 
the specific heat on temperature. In Fig. 2 two typical 
dependencies AC(T) are shown for 6 = (red dashed 



curve) and for 6 = 1/2 (blue solid curve). In addition 
to the differences between plotted curves near the criti- 
cal temperature one can see that in the superconducting 
state at 6 = 1/2 the AC has shallow peak due to the 
interaction of two critical modes. 

Note that the exact solution for the specific heat 
in two-modes critical regime is in a sharp contrast 
with the predictions obtained within the independent 
modes approximation. For comparison in Fig. 1 and 2 
we have plotted corresponding dependencies calculated 
without considering the modes' interaction (black dash- 
dot curves). One can see that these curves differ quali- 
tatively from the ones corresponding to exact solutions. 

IV. CONCLUSION 

Thus we described the critical fluctuations in the 
regime when two interacting modes are strongly fluctuat- 
ing. Both specific heat and magnetoconductivity are re- 
veal special resonant behavior for magnetic flux 6 ~ 1/2, 
which may serve as clear indication of the two-mode crit- 
ical regime for eventual experiments. Note that our re- 
sults obtained for superconducting ring can be general- 
ized directly for superconducting disk^i or the supercon- 
ducting layer with the hole of the round shaped. In these 
systems the magnetic field which exceeds the upper crit- 
ical field H C 2 but is less than the field of the surface 
superconductivity H c % leads to appearance of the super- 
conducting nucleus localized in the ring with the width 
of the order of £o near the superconductor edge. 

The predicted two modes critical fluctuations regime 
produces a relatively strong variation of magnetoconduc- 
tivity and specific heat (of the order of 30%-100% near 
the half flux field). The first detailed measurements of 
the paraconductivity in small superconducting loops have 
been reported already in Ref. [Tl|, and recently the precise 
contactless studies of the magnetic response in individ- 
ual normal metal rings have been performed^. Therefore 
the experimental studies of the critical regime in magne- 
toconductivty seem to be quite feasible. Five years ago 
the sensitive atto joule calorimetr y 18 ! 19 revealed the spe- 
cific heat oscillations in magnetic field for large arrays 
(450 thousands) of noninteracting micrometer-sized su- 
perconducting loops. In such experiments the averaging 
effects may play an important role. Unfortunately the 
statistical characteristics of the studied ensemble of loops 
are unknown, which does not allow comparing their re- 
sults with our calculations. However the recent progress 
in this domain should permit to perform the calorimetric 
measurements on a single superconducting ring in near 
future. 
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APPENDIX 



In the Gaussian regime one can calculate the part of the conductivity which can become singular near the tempera- 
ture of the superconducting transition for < < 1/2, considering only terms with m = and m = 1 in the expression 
©. Then using the fact that in the Gaussian regime for k ^ I the identity (|^ fc (0)| 2 |^(0)| 2 ) = (\i/j k {0)\ 2 ) (\ipi(0)\ 2 ) 

OO 

is true and that J exp (— £/T n ) dt = r n wc obtain 



a ° = T (^r) 



2 t o (IV'o(O)I 4 ) - 0r o (|V'o(0)| 2 ) (» - <t>) (l^(O)l 2 ) 



Tl#0 



(l - 0) 2 n (|^(0)| 4 ) + (l - 0) n (0)| 2 ) <IV^(0)| 2 ) 



Substituting the expressions ([6]) and ([7|) into (|27|) we obtain 



(27) 



1 / e \a (4m£ 2 T 



20 2 t o (frr 



2(1 



7"! 



+ 



(i - 0) n 



E 



(n-0) 



* ^* + (WW™-20)] 



E 



(n-0) 



(28) 



To calculate the principal values of the sums in the expression (|28|) near the critical temperature one can put e* = 
in all terms except those with n = 0,1. The results are 



£* + (6/i?) 2 (l " 20)f [e* + (£o/i?) 2 (l " 20)] ^ [e* + {^/R) 2 n(n - 20)] 



(1- 



E 

(1 + 



(n-0) 



£* + (to/R) 2 n(n - 24,)] 



s* + (Ho/R) 2 (i-24>) (£oAR) 2 U + 20) 



E 

Tl>2 



(n-0) 



(n + 



(e /i?) 2 n(n - 20) (£ /R)*n(n + 20) 



(1 



(1 



and 



e* + (Wi?) 2 (l-20) (eo/i?) 2 (l + 20) 
(n - 0) (1 + 



20 



(£)" 



1 - 120 2 



E [ £ * + (e o /i?)2 n ( n _ 2 0)] 



R 



■ 2 (1 - 40 2 ) 
1 - 120 2 



- — ctg (2tt0) 



40(l-40 2 ) 2 



(G) 

Summarizing wc obtain the final result for the singular part of the component <7q in the Gaussian regime: 

y (1-0) 2 

(1 



--ctg(27r<; 



aP(<f>,e* 



£o\ 4 2 



&Y 



4tts 2 \Rj (£*) 



Co 



Co 



4tts 2 Vi?7 (£*) 3 ' 27ts 2 \Rj [e* + (^ o /i?)2(i_20)] 3 

/ 1 - 120 2 71 



(1 



i?7 [e* + (£o/i?) 2 (l - 20)] (1 + 20) V40(l- 40 2 ) 2 



Ctg(27T0) 



(1 



47TS 2 \ RJ [e* + (e o /i?) 2 (l-20)]^ 



£o\ (1 + 



i?y e* (1 + 20) V 4 0(!- 4 2 ) 2 



1 - 120 2 7T 



Ctg(27T(J 



(29) 



(30) 



(31) 
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For magnetic flux values which are far from the point 0=1/2 the analogous singular part of the conductivity in 
the critical regime can be calculated under the assumption that the fluctuations of only the mode with m = are 
critical while the fluctuations of other spatial modes are Gaussian (the lowest mode approximation). The procedure 
of calculation is similar to the one for the Gaussian regime. The expression takes the form 



-(—)' 

T KmRJ 



WhM0)| 2 / \Mt)\ 2 dt)-cf>( / |^oW| 2 *)^(n-<A)(|V„(0)| 2 ) 



(l - 0) 2 Tl <|Vi (0)| 4 ) + (1 - 0) n <|Vi (0)| 2 ) £(n-0) <|V^(0)| 2 ) 



(32) 



The first term in the square brackets in the expression ([32]) represents the central result of Ref. ll_3j. To calculate the 
second term one should take into account that in the critical regime 



\Mt)\ 2 dt 



4Rs Gif 0) e* ' 



(33) 



The last sum in the square brackets in the expression (f3"2"j) is similar to the one in the expression (f2"T|) except the term 
with n = 0, which should be calculated in the critical regime. The expression for the corresponding correlator in the 
explicit form is given in Ref. [ill Summarizing we obtain the resulting expression for the conductivity in the critical 
regime for magnetic flux values which are far from = 1/2: 



2tts 2 



7T 3 /2 S 2 I R 



* 1 -ml 



Gi 



R I Gi 



1 -ml 



(0) 



(0) ' 

(1- 



2ns 2 V R 



(1 



(1 



+ (do/Rf (l - 

/ 1 - 12(f) 2 



k* + (eo/i?) 2 (l-20)] (1 + 20) ' V40(l-40 2 ) 2 Ctg(2?r</,) 



(34) 



6) 



2tt 3 / 2 s 2 \ Rj Gi {0) 



0(1-0) 



e* + (io/R) 2 (1 - 20) 



Thus the expressions pip and (|34[) give the singular part of the conductivity ao in the Gaussian and critical (in 
the lowest mode approximation) regimes respectively. 
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